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The spinon continues to be an elusive elementary excitation of frustrated antiferromagnets. To solidify evi- 
dence for its existence, we address the question of what will be the Angle Resolved Photoemission Spectroscopy 
(ARPES) signatures of single crystal samples of Herbertsmithite assuming it is described by the Dirac spin liq- 
uid state. In particular, we show that the electron spectral function will have a linear in energy dependence near 
specific wave vectors and that this dependence is expected even after fluctuations to the mean field values are 
taken into account. Observation of this unique signature in ARPES will provide very strong evidence for the 
existence of spinons in greater than one dimension. 



I. INTRODUCTION 

Frustrated quantum magnetism is of active cuiTent interest 
because of the possibility of supporting exotic ground states 
and excitations. One of the central issue that is driving the 
field of frustrated magnetism is the search for such exotic ob- 
jects in Nature. A prime example of an exotic excitation is a 
Spinon - a spin- 1/2 fermion that is charge neutral - which is 
also the main object of our study. The existence of a spin- 1/2 
fermion particle that can be excited out of a ground state of or- 
dinary spins would indeed be remarkable; thus, it is important 
and appropriate to flesh out what would consitute as concrete 
evidence of the existence of the spinon. 

In one dimensional systems, there is substantial evidence 
for the existence of spin 1/2 solitons', though these are dif- 
ferent types of excitations than the spin 1/2 spinon proposed 
in frustrated magnets^. One remarkable source of this evi- 
dence is from ARPES experiments performed on SrCu02^^, 
a tool which has been successfully applied to a number of 
Mott insulators"^. Even though ARPES experiments require 
single crystals in order to have good momentum resolution, 
they are feasible with smaller crystals than what is required 
for neutron diffraction studies. ARPES may also thus provide 
an additional avenue to falsify any of the various ground state 
proposals in the frustrated antiferromagnets. 

One exciting candidate material for a future ARPES experi- 
ment is Herbertsmithite (ZnCu3(OH)6Cl2) which has got a lot 
of attention in recent years. The absence of magnetic order- 
down to 50 mK despite a Curie-Weiss temperature extracted 
from susceptibility measurements of around 300 K makes it 
one of the promising candidates for a quantum spin liquid 
in two dimensions. Herbertsmithite is a layered Mott insula- 
tor with the Copper atoms arranged in a Kagome lattice each 
carrying a spin 1/2 moment much like in undoped cuprates. 
As is well known the Kagome lattice is the most frustrated 
lattice having the least number of neighbours among two- 
dimensional lattices with triangular plaquettes which provide 
geometrical frustration. In multiple magnetic susceptibilty 
and heat capacity experiments^"^, the strong magnetic field 
dependence of the linear-in-(low)temperature specific heat 
suggests that the heat carriers are spin- 1/2 fermions thereby 



making a possible case for Herbertsmithite being a gapless 
spin liquid. Also, neutron scattering measurements on pow- 
der samples have found no evidence of a spin gap''. Perhaps 
the strongest evidence for spinons, however, comes from a 
recent neutron scattering study^ on single crystals'"'" which 
found evidence for a spinon continuum but little evidence of 
a threshold for spinon production in momentum space. If 
ARPES could be applied to this material then it could provide 
complimentary evidence for spinons and solidify the case for 
their existence. 

There is also substantial evidence for spinons in Kagome 
antiferromagnets from theoretical calculations. Recently, 
Kagome antiferromagnets have been shown to have a Z2 spin 
liquid ground state in large scale DMRG calculations with a 
gap to all excitations'^. This strongly suggests that the exci- 
tations are spinons. However, these spinons may not be the 
spinons observed in the above experiments since they have 
not seen any evidence for a gap. Likely, this implies that 
Herbertsmithite is far enough away from the ideal nearest- 
neighbor Kagome Heisenberg AFM with its putative gapped 
ground state through the presence of perturbations such as 
those due to Dzyaloshinskii-Moriya (DM) interactions and 
impurities^"^. Two natural candidates for a description of the 
phenomenology of Herbertsmithie are then a quantum critical 
point on the verge of magnetic ordering arising from the DM 
interactions'^ and the gapless Dirac spin liquid(DSL) stat o'"'''^ 
that is strongly supported by VMC studies''' '*. 

In this paper, we study the effect of Dirac spinons on 
the signatures of a Angle-Resolved Photoemission Spec- 
troscopy(ARPES) on Herbertsmite. Even though the prop- 
erties of Herbertsmithite's ground state are still unknown, we 
will focus on the DSL or its generalization in the form of Al- 
gebraic Spin Liquids(ASL) both because it is a candidate de- 
scription for Herbertsmithite and by itself is a paradigmatic 
state of matter with generic critical correlations. For example, 
the ASL state has also been considered in the context of un- 
derdoped cupratesii. Our main finding is that Dirac spinons 
lead to a linear energy dependence in the ARPES spectrum 
for low energies even though one would expect a sub-linear 
dependence from the anomalous exponent of the spinon prop- 
agator that leads to the aforesaid critical correlations. 
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To expose this experimentally falsifiable linear low-energy 
behaviour of the ARPES spectrum, the paper is organised as 
follows: In section I, we describe briefly the model Hamilto- 
nian and its description in the slave boson picture. In Section 
II, we look at the mean-field theory of the Dirac state and cal- 
culate the ARPES spectral function for this state. This sec- 
tion may be of particular interest to experimentalists since it 
contains our key prediction regarding the ARPES signature of 
the Dirac spinons. In Section III, we go beyond mean-field 
and see how fluctuations of the mean-field do not modify the 
energy dependence of the mean-field ARPES spectral func- 
tion. In our study, we have focused on small fluctuations of the 
phases of the mean-field and thus our study is incomplete with 
regard to large (0[27r]) phase fluctuations. Technically speak- 
ing, the gauge field that models these beyond mean-field phase 
fluctuations is non-compact in our calculation. We conclude 
in Section IV by discussing our results - especially pointing 
out the fact that our prediction might be unique to the Dirac 
state - and possible other experiments to pin down the putative 
spin-charge seperated excitations of the spin liquid. 

While composing this manuscript, a recent similar study 
with similar conclusions has come to our attention"^. These 
authors also consider a Fermi surface state relevant for organic 
spin liquid candidate materials that is not discussed here. 



II. MODEL 



where bl is called a "holon" for it creates positive charge and 
/pct is called a "spinon" and carry the spin index. Here, b 
satisfies a bosonic commutation relations while / satisfies a 
fermionic anti-commutation relations and their product satis- 
fies the fermionic anti-commutation relations expected for the 
electron. Inserting these expressions into ([T), the t-J Hamil- 
tonian in the slave-boson language becomes 

{r,r'),{T (r,r') 

- (Me/2) J2 flfr- - ^l^r + 1 ^ (4) 

where S — a' fl^iy,^' fa' ■ The equivalence of this repre- 
sentation of S and the earlier one in terms of c can be checked 
by comparing their matrix elements in the constrained Hilbert 
space defined by Eq. (|2]l. In addition to the new form of 
the Hamiltonian, the constraint now becomes an equality ir- 
respective of the filling 

fljr.'^ + ^l^r = 1 (5) 

at each r. It is now in a form where we can bring to bear the 
machinery of Lagrange multipliers. 



Consider the t-J model on the Kagome lattice. As is usual, 
the underlying model is the Hubbard model in the large-[/ 
limit. The Hubbard model with a single site U term is a good 
starting point for the spin physics of the Mott insulator Her- 
bertsmithite similar to undoped cuprates. The t-J Hamilto- 
nian is 

H — —t ^ ^ ct ^Cf'^o- + J ^ ^ Sp ■ — fie ^ ^ 

(1) 

where the sum is over nearest-neighbour pairs of sites on the 
kagome lattice shown in Fig. [l] tr is the spin index and S = 
J2cr a' cJ-^cr.cr'Ccr' with T the vector of Pauli matrices. This 
model is understood within the constrained Hilbert space of 
one or less electron per site which can be written as 

J2 4,<.Cr,. < 1 (2) 

(7 

for any r. To handle this constraint mathematically, the slave- 
boson approach is often used. This approach consists of 
increasing the Hilbert space and, in the process, converting 
the inequality constraint to an equality constraint. Also, the 
slave-boson approach is naturally tailored to accomodate spin 
liquid ground state ansatzes [Reference: Baskaran, Zhou and 
Anderson]. 

The slave-boson approach starts with "attaching" a boson 
to an fermion to compose an electron in the following way 

Cra = bifpc (3) 



m. MEAN-FIELD THEORY 

As mentioned in the Introduction, we choose Dirac Spin 
Liquid as our mean field ansatz. The quartic terms of the 
Hamiltonian are decoupled using an auxiliary field and writ- 
ten as follows : 

H = -tXO X! ^rr'blbp - JXO Y Sf?f'/i^/r',£T 
- E (/^Z E f^M- + f^f^l^r^ 

where xo, A*/ and /if, are mean field parameters and sp fi — 
±1 according to the phases defined in Fig. [T] The self- 
consistency condition of the mean-field theory is then given 
by 

\{fljr',a)\ = xoforn.n.r,r' (7) 

(/i,/rv> = (8) 

(biMa) = X (9) 

where the expectation is evaluated in the ground state and jji, 
and jif are chosen to impose the latter two constraints and 
replace /i and the constraint of Eq. (|5]l . The DSL mean- 
field ansatz corresponds effectively to a background magnetic 
flux threaded through the Kagome lattice as shown in the Fig. 
[T] It is chosen so that 7r-flux pierces the hexagons and zero 
flux through the triangular plaquettes. As we see in Fig. [T] 
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photoemission remains small. It turns out, to achieve this at 
zero temperature, any negative value of /if, will suffice. How- 
ever, a specific value will be chosen at finite temperatures in 
the experimentally relevant regime. We will therefore leave 
lib as a finite negative valued parameter in our theory. This 
produces the spectrum for the holons on the right side of Fig. 



FIG. 1: In the left figure is illustrated the choice of gauge for the 
phases for the DSL mean field sppi ; thick lines correspond to a neg- 
ative sign while thin lines correspond to a positive sign. This choice 
(same as Ref. 22) leads to a doubling of the unit-cell. Though the 
gauge-invariant fluxes (+1 through the triangles, —1 through the 
hexagons) respect the underlying symmetry of the kagome lattice 
as they ought to in a spin liquid. In the right figure, we show the 
Brillouin zones - the hexagon is the original Brillouin zone and the 
dashed rectangle is halved Brillouin zone due to the doubling of the 
unit-cell. Q and —Q are the momenta at which the spinon chemical 
potential lies at half-filling. The band structure leads to Dirac cones 
at these momenta as shown in Fig. (2] 




FIG. 2: In this figure, we show the mean field band structure for 
the Spinous and Holons. On the left we show the lowest four bands 
(there are two flat bands at higher energy which are not shown in 
the figure). Of these four bands, at half-filling, the lowest three fill 
out thus leading to Dirac cones for the dispersion of the low energy 
spinous. On the right we show the lowest band for the holon which 
is unoccupied as there are no holons at half-filling. 



the flux still respects the symmetry of the lattice even though 
the particular choice of gauge or phases for this ansatz has a 
doubled unit cell. 

Let us now turn to the spectrum of this Hamiltonian. We no- 
tice that the resulting quadratic Hamiltonian of Eq. |6]consists 
of non-interacting holons and non-interaction spinons hop- 
ping in the presence of the same background magnetic flux. 
The only difference that is the energy scale associated with 
holons is t while for spinons it is J. With this in mind, we di- 
agonalize the Hamiltonian following the notation of Fig. 1(a) 
in Ref. ifTsIl where Bravais lattice vectors are chosen to be 
ai = X and a2 = (1/2)^ + {\/3/2)y giving the Brilliouin 
zone shown in Fig. [T| When this quadratic Hamiltonian is 
solved we get the band structure shown in Fig. |2]with Dirac- 
like nodes as the fermi surface. 

Finally, the question of the holon chemical potential ^ib 
needs some attention. We assume that the holon chemical po- 
tential is chosen so that the doping a; = and the bosons 
remain uncondensed at zero temperature. These assump- 
tions appear quite valid for an ARPES experiment on Herbert- 
smithite so long as the density x of induced holes created by 



A. The Spectral Function 

Let us now turn to the mean field predictions for the ARPES 
spectral function defined by 



A{k,u}) 
-ilm 



-ir-k—iujt 



(10) 



At the mean field level, this becomes 



Al ,A2 

(11) 

where x^^X''^ the holon and spinon bands respectively and 
Ai, A2 are band indices. The above expression has a clear in- 
terpretation: all the possible spinon-holon combinations with 
the right energy-momentum values combine to give the requi- 
site hole at wavevector k and energy uj in the ARPES experi- 
ment. The factor of two with the holon or t term in the above 
equation comes from spin counting but we need not overly 
worry about such factors. Ultimately, we are concerned with 
the u or k dependence of the spectral function, and the func- 
tional form of this dependence is not going to be affected by 
such factors of two. 

We are interested in the low energy behavior of the sys- 
tem in general and the spectral function in particular. We can 
easily see for the mean-field bands, most of the bands won't 
play any role in the low energy behavior It is only the high- 
est filled spinon band and the lowest empty holon band that 
will contribute to the cj — >■ Umit. For these bands, the low 
uj contribution will come from near specific wavevectors that 
connects the "Dirac" points of spinon bands to the minima of 
the holon bands. Near these specific wavevectors, the spinon 
and holon spectrum are linear and quadratic in \q\. Thus, for 
smaller and smaller u, the linear term starts dominating in 
the integral and we get Am^an oc oj for w — > 0. We show 
in Fig. [3] a numerical confirmation of the mean-field result. 
In addition, a back-of-envelope estimate tell us that linearity 



x'iQff 



where Q f and Qb are the 



is expected up to about x'iQb) 
momenta at which the Dirac point and the holon minimum sit 
respectively. This linear behavior in energy is — if the mean- 
field theory were qualitatively correct — a key signature easily 
extracted from ARPES data. 
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A{oj) 




near the spinon chemical potential) - a linear combination 
of tpx, details of which can be found in Appendix A of Ref. 
- and the holons are described by a non-relativistic complex 
scalar field (p (due to the quadratic spectrum near the holon 
chemical potential). The Lagrangian for this field theory is 



L 



free 



drdt 



*(r,t)7^9^*(r,t)+ 



'^holon 



(14) 



FIG. 3: In this figure, we show a numerically calculated represen- 
tative spectrum for ARPES at the mean-field level. The numerical 
calculation confirms the linear behavior at low energies (highlighted 
using the rectangle) as argued for in the text (Sec III). Our choice of 
parameters were J = 0.5 and t = 1. 



B. Continuum limit and scaling analysis 

To understand fluctuations around mean field theory, a con- 
tinuum limit approximation is essential. To this end, it is use- 
ful to carry this out first directly at the mean field level. Such 
an analysis will thus facilitate the scaling discussion when we 
go beyond mean-field theory. We define the continuum fields 
through 

11 

fiiA.o = Y.'l^m,xAR)UTy'^7-'^^^ (12) 



where R refers to unit-cell index and d refers to the index of 
the bases within the unit-cell. The sum is understood to be 
restricted to the low-energy modes, n refers to the index char- 
acterizing the (three) holon minima in the holon band struc- 
ture and are the momenta in the Brillouin zone where the 
minima sit respectively, m refers to the index characterizing 
the (two) Dirac points in the spinon band structure and 
are the momenta in the Brillouin zone where the Dirac points 
sit respectively, a is the spin index, and A is the "Dirac" in- 
dex characterizing the branches of the Dirac cone. and 
C/™^ are the eigenmodes of the bandstructure corresponding 
to holon minima and Dirac points respectively. 

Thus, using the definition of Eq. [TO] we can get for the 
continuum spectral function 



A{QJ - Ql + q,uj) 



A,cr 



(13) 



where 



A,<7 



are contributions due to the eigen modes 
and U™^ ^Ix'a 1^ independent of momentum and energy, 
and contains all the momentum and energy de- 

pendence which is what we will focus on in the following. 

With the definitions in Eq. [121 the low energy physics (of 
the quadratized Hamiltonian Eq. |6|l is described by an ef- 
fective 2+\ continuum field theory where the spinons are de- 
scribed by a Dirac fermionic field (due to the linear spectrum 



The mean field fixed point is then determined by how these 
fields scale when space-time is scaled (and from this field scal- 
ings, we can write a scaling form for the Spectral function). 
When we scale as r — > 6r; t — > ht, it can be shown that scal- 
ing the fields as — > 6^^^*; b^^(t> leaves the Lagrangian 
invariant. Using these field and space-time scalings and the 
continuum limit of Eq. [TOl 



A-mean {k, io) CX 



Im / dr 



dt{(l){r, t)-^{r, i)*(0, 0)0* (0, 0))e-*'- 



ir-k—iujt 



(15) 



we can show that spectral function scales as 
A„^ean{b^^^^,b^^k) = 6"Mmeo«(w,fc) which is Consistent 
with Amean « w as wc obtained from more elementary 
considerations. 



IV. BEYOND MEAN-FIELD THEORY 

The t-J Hamiltonian remains invariant under a local gauge 
transformation made to the spinons and holons 



bij fi.a 



^ fi.a 



(16) 



Thus, when we include the fluctuations beyond the mean 
field level, the resulting continuum Lagrangian should remain 
invariant under this transformation. As is well know, this 
is achieved through replacing the derivative by a covariant 
derivative — >■ I?^ — — ieA^ where e is the electric 
charge and A^ is a Maxwell vector potential satisfying the 
usual transformation law: A^{r,t) — > A^{r,t) ~ d^9{r,t) 
when *(f,i) ^ e*^(^"''*)«'(r,t) and (/)(f,t) ^ e^^(*''*)0(f, i). 
Thus, the beyond mean-field Lagrangian for the spinons and 
holons is 



drdt 



^*{r,t){iDt + {Dl + Dl)l2m- ^x)4>{r,t) 



(17) 



where F^^, = d^A^ - d^A^^. 

Now, Dirac fermions coupled to an U{\) gauge field in 
2 + 1 dimensions is a strongly-coupled problem. To make 
progress, the gauge field contributions are typically handled 
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by a large-N perturbation serie a^^'^^" — arising from N flavors 
of the Dirac fermions. In this case, for large N, 1/N pro- 
vides a formal small parameter provided the coupling scales 
with as ^ When the pertubation theory is carried 
out, one finds that the spinon propagator acquires an anoma- 
lous dimension because of the gauge field coupling e (See e.g. 
Ref. lH^'s Appendix B). Formally, when we scale space-time 
as r — )■ hr\t ^ bt, the Dirac field scales as ^I^ — > 
where S is the anomalous dimension at the fixed point. This 
anomalous dimension is the manifestation of the fact that the 
gauge-field contributions make the theory flow to a new fixed 
point (with different scaling) away from the mean-field fixed 
point. The general intuition is that the strength of fluctuations 
will only increase when TV is decreased and thus this new fixed 
point (called Algebraic fixed point) may dominate a large por- 
tion of the low energy sector for the physical case of = 4. 

Our aim is to understand the effect of the beyond mean- 
field fluctuations on the holons, a subject not visited in previ- 
ous studies due mostly to a focus on the Heisenberg model. 
In an earlier worbi?., the electron spectral function for the 
Dirac Spin Liquid was considered in the context of under- 
doped cuprates, but this study differs significantly from ours 
in that they assumed that the bosons were condensed as they 
were motivated by pseudogap phenomena at finite doping. 
To reiterate our main assumption, we will suppose that the 
bosons are uncondensed, i.e. /i < 0, an assumption that seems 
reasonable if a successful ARPES experiment is to be per- 
formed. In the ideal case at half-filling, under the exact single 
occupancy constraint, there are no bosons around. When an 
ARPES photon ejects an electron, a boson is created before it 
gets destroyed by the replenishing of the electron through an 
electronic bath connected to the sample. This boson is then 
uncondensed at the temperature of the experiment because 
they are very dilute. If the electrons fail to be replenished, 
however, the bosons density could grow to a level where they 
condense. The resulting charged surface, though, would then 
affect the trajectory of the ejected electron and the ARPES 
experiment would fail to obtain meaningful data. Hence, a 
successful ARPES experiment would probe the undoped in- 
sulator only when the bosons remain uncondensed. 

In the following Section IV. A, we discuss in detail the ef- 
fect of gauge field on the bosonic field. We find that the gauge- 
field renormalized bosonic propagator does not acquire an 
anomalous dimension unlike the spinon propagator and still 
scales as a free boson, b^^cj) upon r ^ br;t ^ bt. 
This is expected due to the finite energy cost of the bosonic 
excitations when they remain uncondensed. Furthermore in 
Section IV B, we show that the electron propagator made of 
the spinon and boson also retains its mean-field scaling upon 
taking lowest order gauge-field corrections. This is intriguing 
when contrasted with the anomalous scaling of the spinon. 



the [/(I) gauge field: 

if^ = / d^rdtcp*if,t){iDt + 



^2 

^ - fib)(l>{r,t) (18) 



2mb 

where = 9^ — ieAf^ {r, t) / sqrtN as stated earlier though 
now extended to arbitrary N , the bare mass parameter mi, 
characterises the curvature of the quadratic dispersion at the 
bottom of the holon band and /i^ is the bare holon chemical 
potential. With the assumption of uncondensed holons, fif, is 
strictly negative, i.e. less than the bottom of holon band. 

To understand the holon propagator in the large-N limit, it 
turns out we need to first understand the photon propagator. 
The bare photon propagator receives a finite renormalization 
due to both spinon bubbles and holon bubbles even at leading 
order in 1 /N. However, the holon bubble does not affect the 
renormalization of the most singular contribution arising from 
the fermionic spinon vacuum bubble^"*. Thus, the 0(1/N) cor- 
rected photon propagator is (see Eqn B5 of Ref. ll26tl ) 



vwwv 



N\q 



(19) 



where the ellipsis represents subdominant corrections that are 
smaller than the term shown as q ^ and in the following we 
will use the Lorentz's gauge, ^ = 1. This Coulomb-like form 
is then the natural interaction that the holons feel. 

For the perturbation theory we follow standard texts^' and 
define the original (bare) Lagrangian with bare parameters in 
terms of the physical parameters (to be determined by the ex- 
periment). We start by defining a renormalized bosonic field 
in terms of the bare field 



(po = Z^^^(f) ■ ■ ■ Z = field rescaling of bosons 
Putting this in to Eqn. [18] the Lagrangian looks like 



(20) 



Jf^ = / d^rdt 



L>2 + D2 

0*(f, t){iDt + ^- ^ - fi)c^{r, t) 

Zm 

2dm 

+ dzrir,t)iDt(bif,t) + 6^^*ir,t)(bif,t) (21) 



where 5z = Z-1, 5,„ = - and 5^ = -Z^ii, + ^ 

with /i a renormalized version of /ib not to be confused with 
the chemical potential of the electrons. These are the countert- 
erms. Since the holons do not renormalize the photon field, 
it's propagator stays the same as in Eqn. [19] 

We rewrite the renormalized Lagrangian as a free part 



free 



d^rdt(j)*{r,t){idt 



dl 



2m 



y)<t>{r:t) (22) 



a counter term part 



A. Renormalized Perturbation Tiieory for Holons 

Let us now focus on the holon part of the Lagrangian — a 
non-relativistic complex bosonic field — that couples only to 



= / d^rdt 



^{r, t)^^-^<l>{r, t) + S^r{r, t)cp{r, t) 

ZOm. 



(23) 
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and an interaction part 



Z^*{f,t)Atir,t)(bir,t) 



^^^^^^^^^ .dMr,t) + dyAy{r,ty ^^^^^^ 



nib 

( Al{r,t)+Al{r,t)\ 



(24) 



Therefore, the Feynman rules are given by: the boson propa- 
gator Hne 



1 



w- (|fc|2/2m+ 1^1) +k 



(25) 



the counter terms by 



and the interaction vertices by 

• Boson-Gauge Field's "Time" component vertex 



(26) 



= Z (27) 
• Boson-Gauge Field's "Space" component vertices 



Zki/nib 



(28) 



where i refers to the Gauge field's "space" component 
we are looking, i.e. if cf) interacts with A^, then the ver- 
tex will contribute k^, etc. The are actually two such 
vertices of this type, one where ki goes with the incom- 
ing boson's ((p) momentum and one where ki goes the 
outgoing boson's (ip*) momentum. We can collected all 
these vertices together by defining V^'' (fc) to be Eq. (|27] | 
for ^ = t and Eq. (l28l l for fi = x or ji = y. 



(29) 



• Boson-Gauge Field's quartic vertex 



-Z/2mb 



Note: energy-momentum is conserved at the interaction ver- 
tices. 

Now, we look at the 2-point free vertex function Fg (fc, w) — 
uj — (|fcp/2TO + |/i|) and impose our renormalization condi- 
tions on the renormalized 2-point function 

r2(fc = (o,o),L.-Af) = M = -l/^l 

drVdw(fc = (0,0),cj = M) = 1 
drVrffc^(fc = (0,0),tJ = M) = -l/2m (30) 



.9) 




+ 



+ 

+ 
+ 



FIG. 4: Perturbation Series for Bosons to first order in 0[1/A'^] 



In other words, we chose the time-like scale (fc,a;) = 
((0, 0), A/) as the renormalization scale. 

With the above setup, let us now turn to our main interest, 
the holon's anomalous dimension which is governed by the 
holon field rescaling Z^'. We will consider this to order 1/iV 
through the diagrams of the propagator shown in Fig. |4] If we 
view the self energy as the contribution to F'^ of r'^(fc. lS) — 
rg(fc, lS) — E(fc, lS) + counter terms, then from the diagrams 
in the figure, we see that the second renormalization condition 
becomes 



1^— ((0,0),Af) 
auj 



(31) 



So, to understand Z we need only consider dTi/duj. Since 
the contributions of the last two counterterm diagrams and the 
diagram with the "q = 0" bubble in Fig. |4]are independent of 
ui, these terms are then unimportant here. Hence we obtain 



dS(fc, uj) 



d f d'^qdqo 16 
doj 



duj dw J (27r)3 N\q\ 

l + ^{2k^-q^.f + ^{2ky~qyf 

b b 

[iuj-qo)-{{k-q)y2m+\fA)+ie)] 
Taking the derivative, we then obtain 



d^,,^ „x 16 /■ d^q 1 



(32) 



(33) 



1 



r{qy 



[qo-ie- M +\q\'^/2m+\^i\ 



To understand this integral, we need the following identity, 
whose origin lies in the branch cut singularity: 



lim 

= 2 



dx 



dy 

(iib + y) - aY^y'^ +2by'' 



(34) 



& > 



Using this, we then obtain 

— ((0,0),M) = A / dq^dqy dy (35) 

J-oo Jo 

''(z(y+|<7l)-M+|glV2m+|M|)2 
1 



Finally, to understand the behavior of the above integral in the 
infra-red, we introduce an UV cut-off A for the integral and 
express everything with the scale of energy in units of this UV 
cut-off A. 

A'^^A' T''^'"T'™" A' ^ ^ 
In terms of these rescaled variables, we then obtain 

— {{0,0), M)^AJ ^ dq^dqy dy (37) 

b b 

(z(y+|||)-M+||2/2m+|/~i|)2 
1 



y2 + 2|q|y 



Consider the infra-red limit lim Af — 0, i.e. when the renor- 
malization scale is in the infra-red. The self-energy integral 
simplifies to 

TV^ /• 1 /• 1 

— {{0,0), M) = AJ ^ dq^dqy dy (38) 



(»(y + k1) + kT/2m + IAI)2 
1 



y2 + 2||y 



The only singularity in the integrand comes from ?/ — ^ 0, 



\q\ — because of yy^ + '2,\q\y. The presence of \fl\ au- 
tomatically regulates the other denominator Let's say we do 
the y integration first. Since, y will dominate over y^ when 



y — >■ 0, therefore y + 2|g|y w y 2|<7|y. As an integral over 
y, this is an integrable singularity. Similarly as an integral 



over (radial) \q\, the y 2\q\y is an integrable singularity. 

Thus, in the infra-red limit M — > 0, d'E/duj is finite 
and hence can't have a logarithmic contribution(log A/ — 
logM/A). Through similar arguments, we can show that 
T,{uj) is finite too. Also, the integral coming from the third 
renormalization condition in Eqn. |30] which determined Sm 
will be finite as well. We mention this because in integral for 
dS /du!, we have the term Z'^/ml in the numerator and mi, de- 
pends on ni and 6,n. But the finiteness of '6,n' integral guar- 
antees us that there wouldn't be any logarithmic contribution 




FIG. 5: The terms in the perturbation series for the Hole propagator 
uptoO[l/7V] 



to dTi/duj because of the presence of the jm^ term. Thus, 
our earlier conclusion - that dY,/duj has no logarithmic diver- 
gence - remains valid. In absence of logarithmic divergences 
the bosonic field won't acquire an anomalous dimension^, 
and thus Z = 1. 



B. Large-N perturbative calculations for an electron 

In the previous section we showed that large-N corrections 
to the holons are finite and thus the holon propagator is 
just a renormalized form of the mean-field holon propagator 
As has been shown before by others, the spinon propaga- 
tor does get a logarithmic correction. In this section, we work 
out the consequences of these two statements on to the elec- 
tron propagator made of the a holon and a spinon that would 
be measured in an ARPES experiment. 

The Feynman diagrams for the electron propagator to 
0[1/A^] are shown in Fig. |5] We are interested in finding 
out if the hole propagrator recieves any logarithmic correc- 
tions due to the gauge fluctuations due to divergent diagrams. 
Let us first look at Fig. |5] a) which is the continuum limit 
of the mean-field electron propagator Amean {q, i^) studied in 
Sec. Ill and shown to be linear in co through an analysis of 
mean-field band structure of spinons and holons and through 
a mean-field scaling argument. For our low-energy continuum 
theory (we use the same conventions for the spinon propaga- 
tors as in Ref. 26), this diagram is equal to 

Aa{q,uj)^ j duJidqiS{qi,uJi)G{q- qi,uj - uJi) (39) 

where is matrix spectral function that enters Eq. [T3]and is 
related to the measured spectral function through matrix ele- 
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ments, 



and 



(40) 



(41) 



is the boson propagator entering Eq 
frequency integral we arrive at 



Evaluating the 



Aa{q,uj) = TT dq 



1 ■ 91 



(42) 



In the infrared limit at the Dirac point w — > 0, (f = (0, 0), 
the integrand in Eq. |42] has a finite denominator, thus there 
is no divergence. This is as we would have expected of the 
mean-field hole propagator and the leading term (in lo) of the 
(finite) integral is linear in a; (as concluded from other argu- 
ments in Sec III). Since in the previous subsection we showed 
that the holon recieves only finite corrections at 1/7V, the di- 
agrams Fig. |5]d), e) and f) are also finite as can be shown in 
a similar manner as the evaluation of the mean-field diagram 
above. Hence, diagrams a), d), e) and f) are all finite diagrams 
as far as divergences are considered and are consistent with 
the mean field scaling. 

We proceed to look at the diagram Fig. |5]b) and c) where a 
photon is exchanged between the holon and the spinon. This 
diagram equals to 



Ab(g, w) — j d(jJidqiduj2dq2 



S((fi,IUi)S((fi - (f2,Wi - W2) 

X G(g - qi,u} - a;i)G'((f- qi + q2,uj - uji + UJ2) 
rD^,iu:2,q2)[V''iq- qi) + V^q - gl + 92)] (43) 



where Df^^, {q, to) is the photon propagator of Eq. ( fT9] l and 
V'^{k) is defined in Eqs. ( |27T i and (|28T l. Doing the oji-integral 
as a contour integral by closing on lower-half plane we arrive 
at 



167r 
If 



dqiduJ2dq2 



«7 - 7 • Qi 



LO + i\qi\ 



2m 



IMI 



7°(~^|gi| - ^2) + 7 ■ (gi ~ '?2) 

{u^ + ^\q^\+U2) + ^e- \^-%+f^\' -\^,\ 
7^[V^(g - gi) + V^^iq- gi + 92)] 



(44) 



and using Eq. |35]to do the uj2 integral this further reduces to 



327r 

IT 



dqidq2 



17 



7 • gi 



UJ + i\qi\ - 



|g-giP 

2m 



X /'"'dy ^°HI'?i|-*(l92|+2/))+7-(gi-g2) 



{uj + i\cii\ + i{\q2\ + y)) 



X Yi^iq- gi) + V'iq- qi + 9*2)] 



\q-qi+q2\^ _ 1,, 
2m l/^l 

1 



^y^ + 2y\q2\ 



(45) 



Looking at w — > 0, q = (0,0) again (the infrared limit at 
the Dirac point), we see that the only singular piece in the 



integrand is 



which (as argued in Sec III b) is an 



^/y^+2v\q2\ 

integrable singularity. Thus, the diagrams Fig. |5]b) and c) 
give only finite corrections to the hole propagator. 

We finally come to the diagram Fig. |5]g), where the loop 
on the spinon propagator has been shown to give a logarithmic 
correction to the spinon propagator We want to know what is 
the effect of this logarithmic correction on the electron propa- 
gator. This diagram equals 



37r2iV 



duJidqi 



(70a;i + 7 • gi)log(^w^ + |gi|VA) 



|gil 



(46) 



where the logarithm is due to photon-spinon loop (see Eq. 
C9 of Ref. |26i). We can express log(Vwf + |gi|VA) as 
(l/2)(log((tJi + zgi/A) + log((wi - ^gi)/A)). As a func- 
tion of (complex) uji, we show the branch cut structure due to 
the logarithms in Fig. |6l Thus we will do the cji integral as a 
contour integral keeping in mind that for the log((aji +igi/A) 
term, we close the contour on the upper-half plane, while for 
log((aji — iqi/A), we close the contour on the lower-half 



plane. By capturing residues we get that the diagram equals 



A 



A + B 
4 



B = 



C 



Z-kN 
4 



37riV 



log(- 



C . 



where 



rfgi 



(-*7" + 7-gi)iog(-^) 



U) - 



«|gil 



\q-qi 



(*7" + 7-gi)log(2iiM) 



dq- 



'-*|gil 
(7°(^ - 



2m 



- m 

ImI) + 7-gi) 



*|gil 



2m 



|gi| 



A 



(47) 
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In the infrared limit at the Dirac point (i.e. w — > 0, (f = (0, 0)), 
the above simpHfies to 



= A + B + C ... where 



A 



B 



3nN 
4 

3nN 



dqi 
dqi 



H7''+7-'?i)(log(^)-f) 



OJ + i\qi\ - 



2m 



(z7° + 7-gi)(log(^) + f ) 



C = 



3ttN 



dq 



OJ - t\qi\ - 
0/, , _ \qi\'^ 

2m 



91 



(7°(^ 



|/i|) + 7 • (fi) 



X log( ) 



ill 

2m 



191 1 



A 



(48) 



The 7 • (ji and 7 • ifi terms are odd under inversion and hence 
upon integration equal zero. A + B is finite. This can be seen 
by expanding the denominators of A and S in jgij. At low- 
est order for \qi\, the two (singular) logarithmic contributions 
cancel each other out. Higher orders in |gi | don't have singu- 
lar integrands as \qi\'"' log(|(fi|) — > as |(fi| — > for n > 1. 
We can write the logarithm in C as 



-Mr+\Qi\ 



■) + iarctan(- 



191 



We 



log( 

'^^^ — ifj-i 

focus on the imaginary part of C since it is the most singular 
part. We can write this as 



3ttN 



dqi 



7°(^-lf-|Ml) 



(,, I'll' 



X log( 



2m 



-) 



A 70(..-H-y)log( ^-l^l'/+^'"^ ) 
may 



(w - 1^1 - + 2my 



(49) 



and we see that the integrand is finite everywhere on the in- 
tegration domain (with the UV cut-off). Therefore the inte- 
gral cannot have any singular/logarithmic term. Thus we con- 
clude that 0\i/N\ corrections to the hole propagator due to 
the gauge fluctuations are finite. Therefore the 0[1/Af] scal- 
ing of the ARPES spectral function is the same as the mean- 
field scaling. 



V. DISCUSSION 

The key prediction of our work is that the ARPES spec- 
tral function for the putative spin liquid, Herbertsmithite, will 
show a linear low-energy behavior if the ground state has 
Dirac Spinons. This prediction is expected to be unique to 
the Dirac Spin Liquid as the Dirac point-like Fermi "surface" 
is not generic. Thus, other spin liquid ground states might 
generically be expected to have different energy dependences. 
Recently, NMR experiments"^ were performed on single crys- 
tal Herbersmithite and the authors reported that their results 
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H 

a) 





9 

- -( 


► jlgil 


b) 





FIG. 6: We show the branch cuts in a) of log( "'-^l^'-l ) and in b) 
oflog(ii^) 



are not in contradiction with a DSL. This makes the case for 
an ARPES experiment on Herbertsmithite even more exciting 
since, as stressed before, our prediction is a falsifiable predic- 
tion. 

Another unique feature of our prediction is based on the fact 
that DSL is a high symmetry state in momentum space. Thus, 
the linear energy dependence will be most clearly observable 
at very specific momenta (i.e., momenta that connect these 
high symmetry points in the Brillouin Zone coming from the 
Dirac points for spinons and the dispersion minima for the 
holons — Q^)). Observation of these momenta would be 
another strong evidence for a DSL. 

We found that the holons in the Dirac spin liquid get trivial 
renormalizations from the gauge field fluctuations and remain 
as free bosons. On the other hand, as has been shown in ear- 
Her literature, the Dirac spinons do get nontrivial logarithmic 
renormalization from the gauge field fluctuations. This fea- 
ture points to a dramatic spin-charge separation in the Dirac 
spin liquid. 

In this paper, we have discussed how ARPES can provide 
an avenue to prove the existence of spinons in Kagome An- 
tiferromagnet. Existence of spinons would raise two obvi- 
ous questions : 1) How would one detect the corresponding 
charged holons ?, and 2) How would one detect the photons 
corresponding to the U(l) gauge field ? 

Detection of holons can be envisaged perhaps through a 
transport experiment. Since the holons are charged bosons, 
we can expect to probe their existence in transport measure- 
ments. Since the low energy holons occur at particular non- 
zero momenta in the zone, there might be a non-trivial angular 
dependence of conductance in a sample which can be com- 
pared against theoretical predictions. Detection of the "pho- 
ton" is probably a harder task. These photons will of course 
contribute to the specific heat along with spinons and one will 
have to separate the two contributions. Thermal conductivity 
measurements can similarly get contributions from the pho- 
tons. 

In conclusion, frustrated quantum magnets offer an uni- 
verse of interesting objects including the one we focused on in 
this study : Dirac spinons, U(l) photons and non-interacting 
bosons built out of a Dirac spin liquid state. Herbertsmithite 
provides with an ideal candidate material to explore this uni- 
verse and find these exotic objects. Finding the Spinon in Na- 
ture will certainly be a big step in the uncovering of that uni- 
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verse and a successful ARPES experiment on Herbertsmithite and Siddharth Parmeswaran for useful discussions. During 
might go a long way in this regard. the work, SP was supported by NSF grant DMR 1005466. 
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